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In this paper the synthesis of a constant, parameter output, food hark 
control law of constrained structure is set in a multiple objective linear 
quadratic regulator (MOLQR) framework. The use of intuitive objective 
functions such as model-following ability and closed-loop trajectory sen- 
sitivity, allow multiple objective decision making techniques, such as the 
surrogate worth trade-off method, to be applied. For the continuous-time 
deterministic problem with an infinite time horizon, dynamic compensators 
as well as static output feedback controllers can he synthesized using a de- 
scent Anderson- Moore algorithm that is modified to impose linear equality 
constraints on the feedback gains by moving in feasible directions. Results 
uf 'three dill’erent examples are presented, including a unique reformulation 
of the sensitivity reduction problem. 

Introduction 

In synthesizing a control law a designer has two objectives in mind 
maximize performance and minimize the cost of implementation. A simpler 
control law, which is less complicated and less costly to implement than, 
for example, a full state feedback controller, may be preferred by the 
designer. The designer may have a number of structural alternatives in 
mind such as full output, feedback, decentralized control, or low order 
dynamic compensation. In this paper the problem of minimizing the cost 
of implementation is viewed as a problem of selecting a control law among 
a number of alternative control law structures. The trade-offs between the 
performance objective and the cost of implementation can only be examined 
fairly if the free parameters of each structure optimize the performance 
objective. The performance objective itself may be characterized as a set 
of objectives. The objectives considered in this paper are integral quadratic 
objectives of state energy, control energy, model following ability, trajectory 
sensitivity, and subsystem objectives. These objectives are often conflicting, 
therefore a multiple objective approach can be taken. 

Algorithms for constant parameter output feedback controllers have 
been extensively studied [1]-[17]. Few algorithms however, allow the de- 
signer to arbitrarily prespecify the control law structure [6], [7], [11], [14], 
[17] Fewer still consider the multiple objective nature of the problem [13]- 
[15]. Three functional scalarization methods have been applied to MOLQR 
problems - the weighted sum, e-constraint, and goal attainment methods. 
All three methods convert the multiple objective problem into a single ob- 
jective problem A single noninferior solution may then be found by solving 
the single objective problem using any one of a number of nonlinear pro- 
gramming techniques or algorithms that solve the necessary conditions. 
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'Hie objective function to he minimized that will he considered here is 
defined as 
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where Q is a positive semi definite n x n malrix, It is a positive delinite 
in x m matrix, and both are symmetric. In order lor J{G,fu) to he finite, 
the states weighted by the matrix Q mnsl lie rendered asymptotically stable 
by application of a finite amount of control energy. 

I'or any feedback gain matrix G that, stabilizes 1 he dosed loop system 
A — A -- HGC the value of the objective function may be evaluated a.s 

■n<lr n ) = 4l<*» (M) 

- tr { A x ( )Zq } 

where A is the unique positive semi-definite solution of the linear matrix 
equation 

A 1 1\ + KA + Q + {GCfRGC = 0. ( 1-5) 

Thus what may have appeared to be a dynamic optimization problem 
(eqs. ( 1- 1)-( 1-3)) , is now formulated as a static objective function (1-4) 
that is to be minimized with respect to G and K subject to the equality 
constraints (1-5). Finally, the problem may be restated as: Given the 
dynamic system (1-1), initial condition x 0 , and a control law of the form 
u(/) = - Gy(t ), choose G * to minimize the objective function J(G,xq). 

The first-order necessary conditions are found by first forming the 
Lagrangian 

C(G , A, L) - tr{ KxqXq } + tr{[d / K + KA + Q + (GC) T 11GC]L J } (1-6) 

where L £ R nXn is a matrix of Lagrange multipliers. By using gradient 
matrix operations [4], the first-order necessary conditions for G* , I\ m , and 
L* to be optimal are given by 


The weighted sum method is used in this paper. The solution algorithm 
used is a descent Anderson-Moore algorithm modified to move in feasible 
directions. This algorithm was chosen for modification because it has been 
shown by example to be faster than the Davidon-Fletcher-Powell method 
[8], and because proofs of convergence to a local stationary point of the 
cost function of Anderson-Moore type algorithms have been reported on 
[8], [11] (these proofs are contested however in [15].) 

This paper is organized into two parts. In part one, the constrained 
optimal output feedback problem is formulated and a solution algorithm is 
presented. In part two, multiple objective LQR problems are formulated 
and results of examples are summarized. 

1.0 Optimal Controllers of Constrained Structure 

1.1 Full output feedback controllers 

The linear dynamic multiple-input multiple-output plant model is in 
the standard state space form 

x(t) = Ax{t) + Bu(t) ^ ^ 

y(t)=Cx{i) 

where x £ R n is the state vector, u £ R™ is the input vector, and y £ R p is 
the output vector. The outputs are assumed to be the sensor measurements. 
The initial state, or initial condition, is given by x(0) = xq. 
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= 2[AG*CL*C 7 - B r ICL*C r ) = 0 (1-7) 

= A* IS + ISA mT T xqx% = 0 (1-8) 

= A* t K' + K m AS + Q + (G‘C) T AG*C = 0 (1-9) 


where A* - A - QG m C. Equation (1-7) may be rewritten as 

G* = R~ l R 7 A * L * C l ( C L * C T ) “ 1 . (MO) 

The optimal cost is J* (G* , xq) = tr{A "a-oJ’o }• 

1.2 Output feedback controllers of constrained structure 

A full output feedback formulation assumes all the elements of G 
are free to be optimized. Although full output feedback uses only those 
measurements that are available, it does not allow the designer freedom to 
choose feedback paths. Feedback paths may be eliminated by constraining 
elements of G to be zero. In addition, if the set of available measurements 
is too restrictive, dynamic compensation may be necessary to stabilize 
the system. Even if a constrained control law stabilizes the system, the 
addition of a dynamic compensator may improve performance. Dynamic 
compensators may be designed by augmenting the state and output vectors 
and constraining appropriate elements of the output feedback matrix to be 
zero or one. 
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The only difference between the first-order necessary conditions for 
optimality for a full output feedback controller (1-7) -(1-0) and an output 
feedback controller of constrained structure (elements of G constrained to 
be constant) is that only the partials of £ with respect to the free parameters 
of G need to equal zero. The additional necessary condition <7* £ U must, 
of course, be satisfied. 

In order to constrain the structure of an output feedback controller, 
elements of the feedback gain matrix G are constrained to be constant. 
These constraints are a special form of linear equality constraints. Linear 
equality constraints arc useful for other purposes as well as constraining the 
structure of the gain matrix. A robust control law for multiple plants can be 
designed using a multiple objective approach by constraining corresponding 
elements of the feedback matrix to be equal (section 2.4.3). Trade-offs 
between closed-loop trajectory sensitivity and other LQR objectives can 
also be examined by constraining elements of the feedback matrix to be 
equal to each other (section 2.4 .ft). The initial values of some elements 
of the control vector may be specified by realizing that the initial control 
vector is u(0) = -CCrfO), [14]. Calisc and Raman use linear equality 
constraints to impose a degree of modal insensitivity on the closed loop 
system [1G]. 

The optimization problem for the constrained optimal output feedback 
problem can be stated as: (liven the dynamic system ( l- 1 ) , initial condition 
.To, a control law of the form u(l) = -Gy(t), ami the constraint sot 

il --Z (G| / t| (G) - b t = (), / - | in) (M2) 

where l he scalar function h,(G) combination of the elements of 

G, choose G* € R to minimize the objective function J(G,xq). 

1.3 Solution Technique 

To enforce the linear equality constraints (1-12) the following penalty 
function is defined 

m 

7(«) = 5 - k?. (M3) 

1 1 = 1 

The objective function tu be minimized becomes J(G,x o) = J(G,xq) + 
7(G). The Moerder-Calise algorithm [11] may be applied directly to this 
problem formulation. However with the penalty function method the solu- 
tion is never feasible and the stationary points of J are not the stationary 
points of J unless 7(G) = 0. The constraints must be forced to be satisfied, 
and then the resulting feasible solution evaluated. Increasing the weight- 
ing factor i' to a large value may cause problems with convergence of the 
algorithm. 

In order to avoid any problems associated with penalty functions, a de- 
scent Anderson-Moore algorithm is modified to move in feasible directions 
only. A feasible direction is obtained by projecting the unconstrained direc- 
tion onto the linear equality constraints at each iteration of the algorithm 
Define the vector \ as the vector of elements of the matrix G. The linear 
equality constraints (1-12) can also be represented by the matrix equation 

A cX = f>. (1-14) 

If c/jt is the unconstrained direction at step k of the algorithm, then ci* the 
projected feasible direction is [24] 


2) For the constraint A g tJ — A </t/, change the sign of both Ay tJ and A< /*/ 
to ensure that. 

4— - A//*, < 0 (M9) 

<hiii rifiu 

To a first order approximation, this results in a reduction in cost along 
each vector of a basis formed by the unconstrained gains A ?Hj and the 
constraints A//, ; — A Note that, in general, these type ul 'sign change' 
rules only apply when the linear equality constraints /i,(G),z = 1 , . , . , m are 
independent. 

1.3.1 The Algorithm 

The modified descent Anderson- Moore algorithm for computing G* is: 

0) Choose («' (1 G 12 such that A — A — HG(\C is stable. Set i — 0. 

1) Solve (1-9) for A,- and (1-8) for /,». 

2) Evaluate 

A Gi - ir 1 il r K t L,C r {('!.,( ,r - G,, 


3) lls< equations ( I- Id) uml (117) to project. AG, onto the constraint set 
12 ami obtain AG, 

*1) Use rules (M8) and (1-19) to change' signs of elements of AG; if 
necessary, where 

- 2 [RGiCLiC 1 ' - B r KiLiC r ] 


5) Set 


G, + i = G, + a AG;. 


where a £ [0, 1) is chosen to ensure 


•A + l < Ji = tf { /v; XqXq }. 


6) Set i = i + 1 , go to 1). 

Practical methods for determining an initial stabilizing gain Go may be 
found in [9]. The algorithm was programmed in Turbo- Pascal on an IBM 
PC-XT. Step 1 requires the solution of two Lyapunov equations (1-8) and 
(1-9). See Smith [25], or the more efficient algorithm of Bartels-Stewart 
[26]. A method for selecting ot is needed in step 5. The convergence proof 
of Moerder and Calise [11] shows that there exists an a small enough so 
that J t+ i < Ji. Therefore the simplest method is to start with a = 1 
then reduce ct if necessary to ensure a reduction in cost. The algorithm, as 
described, lacks a stopping criteria, however the criteria \dJ jdG{\ < fj or 
Ji — 1 < c 2 may be used to stop the algorithm. 

2.0 Multiob jective LQR Problems 


dk = -[/ - Aj(AcAj) l A c ]<lk- 0-15) 


Define AG to be the unconstrained direction and AG to be the pro- 
jected direction. In this paper the following two constraints and associated 
projections are used: 

1) To constrain the gain element ijij equal to a constant c,-j, the projection 

Ag { j = 0. (1-16) 


2) To constrain g u equal t.o gth the projection is 


A Qij - Agu = ( 


A gij + A<m 
2 ; 


(1-17) 


Two rules were used to change the sign of the elements of AG in order 
to ensure an improvement in cost at. each step of the algorithm: 

1) For unconstrained, change the sign of A g^ to ensure that 


2.1 Problem statement 

The problem is to optimize simultaneously N objective functions. The: 
multiple objective problem, or vector optimization problem, is stated as: 

min {7, (G), J 2 (G), .... Jjv(G)} (2-1) 

where G is a vector (or matrix) of decision variable or optimizing parameter 
values and Q is a constraint set. Each objective function 7,(G) is an integral 
quadratic objective function similar to J(G,x 0). The explicit reference to 
the initial condition xq is omitted for notational convenience. 

If the objective functions are conflicting, which is often the case, the 
global optimums of J;(G), i = 1, . . . , N, are not achieved at the same GV 
Clearly there is no unique solution to the problem. However, a set Q* of 
noninferior solutions (or Pareto optimal solutions), can be defined. 

DEFINITION: A decision G m is said to be a noninferior solution to 
the problem posed by (2-1), if there does not exist another G so that 
Jj (G) < Jj(G*), j — 1,2, ...,7V, with strict equality holding for at 
least one j [18]. 


— Aj ( j<0. (1-18) 

fyij 
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Therefore a noninferior solution is one in which a further optimization 
of any J,(G) is at the expense of at least one of the others, (i.e., all 7, (G)'s 
cannot be reduced simultaneously). The concept of a noninferior surface is 
graphically illustrated by the two objective, single decision variable problem 
in figure 1 where objective functions Jj(<7) and 7 2 (r/) are plotted versus the 
single decision variable//. The objective functions Ji(g) and ^2(y) achieve 


2 



t heir optimums at </" and y ^ respectively. The values of y between g* { and 
tj* are noninforior solutions to the problem because a further decrease in 
either objective leads to an increase in the other. 1 he noninferior surface 
(or t.rade-ofT surface) is plotted in the decision spare in figure 2. 

d ims the problem of finding a solution that is acceptable to the control 
system designer is reduced to finding an acceptable soluLion in Llie notiiii- 
ferior set.. UV rather than the entire feasible set, Q. To solve the overall 
problem then, a trrhmqiie for generating noninferior solutions is needed as 
well as a technique for determining which noninferior solution the decision 
maker (control systems designer) prefers. 

2.2 Generating noninferior solutions 

in this paper the weighted sum method is used. Noninferior solutions 
are generated by solving a single objective problem that is the weighted 
sum of the S objective functions. The problem is staLed as: 


min 

O'tlt 


A' 



0 < < 1 


( 2 - 2 ) 


( ieoim I in ally, for a two objective problem, the solution of (2-2) is 
represented by the closest line to the origin, of slope ~W'j/w\ x that is tangent 
to Hie noninferior surface (figure 2). For any set. of weights satisfying 
V »r, - i the solution to (2-2) is a noninferior solution. However, the 
weighted sum method fails to lind all points on a nonconvex trade-oil 
surface. In figure 2 the noninferior solutions between points A and B cannot 
be found using tin* weighted sum method. The weighted sum method has 
the advantage that if the weighted sum of objective functions maintains the 
same form as a single objective function, an algorithm that applies to the 
single objective problem may be directly applied to the multiple objective 
problem. This will be shown to he the case for MOLQR. problems. 

2.3 Deciding on an acceptable noninferior solution 

Although a nonlinear programming algorithm may be used to find a 
single noninferior solution, the problem of finding a noninferior solution 
that is preferred by the decision maker (the control systems designer) must 
still be addressed. In making a decision it is typically assumed the decision 
maker maximizes a utility function i.e., 

max /7[Ji(G), . . . , Jn(G)] (2-5) 

where / [ j is a monotonic decreasing function of the objective functions 
and Q* is the set of all noninferior solutions. There are many 
interactive techniques that can be used for extracting from the decision 
maker the addit ional information necessary to solve (2-5) [19]. A technique 
which requires a minimum amount of information is the surrogate worth 
trade-oil (SWT) method [18], [19]. The SWT method uses the values of 
the trade-off functions, defined as 


dJi 

dlj 




( 2 - 6 ) 


in questioning the decision maker. The decision maker (DM) is asked a 
irade-ofF question such as: “Given levels of objectives J\(G) X . . . , J^(G) 
ho w w illing are you to trade A*j units of Ji to reduce Jj by one unit?.” The 
DM is questioned about a sequence of noninferior points until a point is 
found such that the decision maker is indifferent to moving in any direction. 


with the weights q XJ > 0 and r*/ > 0 chosen to ensure Q > 0, R > 0, 
and both Q and R symmetric. The weights may be normalized so that 
j ( hj + Yh j r u — I Thus the solution to (2-7) may be interpreted as 
a noninferior solution to the multiple objective problem consisting of Un- 
objectives weighted in (2-8) and (2-9). To arrive at a noninferior solution 
that is acceptable' to the designer the weights art' changed and the problem 
is solved again. However the objective# in (2-8) and (2-9) containing cross 
terms (e g., x t (t)xj(t), i ^ j) in general have no physical meaning to the 
designer. Tims it is dillicilll, if not impossible, for the designer to answer 
the question : “Are you willing to trade A I; units of J, for one unit of J 3 7” 
The problem may be reformulated as follows: Consider objective funct ions 


of the form 



• *o 

f (Xi(t)-X } (t))\ll, 

Jo 

f (2-10) 

Jo 

f ..?«)*, 

to 

f (M0 - »j(0)“d£, 

./i) 

f («,(() + u 1 (t))-dt. (2-11) 


A standard hQH objective function (2-7) is obtained by weighting these 
six terms such t hat Q > 0, R > I), and both Q ami R are symmetric. The 
integrands of the integral equations in (2-10) and (2-11) are positive for 
all /, therefore weighting (he objectives in (2 10) by < f X) > 0 ensures that 
Q > 0 ami weighting the objectives in (2-11) by weights r }} > 0 ensures 
that R > 0 I bis observation motivates the following theorem 

Til FOR KM 2 I 

Given a multiple objective linear quadratic regulator (MOLQR) prob- 
lem that, is defined as a weighted sum ol Objective functions as 


J(G) = 


/ [* T 

Jo 


(t)Qx{t) + u 1 (t)Ru(t)} dt 


= u>i • 


the noninferior surface between the N objectives 7, is convex for a full state 
feedback control law if for every set of weights that, satisfies 

1) 0 < in, < 1,*— 1 N 

2 ) £*, = ' 

the penalty matrices Q ami R satisfy 

3) Q positive semi-defmite ( Q > 0), and symmetric 

4) R positive definite (R > 0) and symmetric. 

PROOF: 

The proof is by contradiction. Assume that the noninferior surface is 
nonconvex. Then there must exist a set of weights such that there are two 
(or more) different G * that minimize 7(G). By assumption, for every set 
of weights satisfying 1) and 2), Q > 0 and R > 0, therefore the minimum 
of 7(G) is unique. Thus there does not exist a set of weights satisfying 1) 
and 2) which results in two (or more) different G* minimizing 7(G). The 
noninferior surface therefore must be convex. D 

The objectives in (2-10) and (2-11) have more meaning to the designer, 
in that they are easier to trade-off, than those of (2-8) and (2-9). 


All iterative methods assume the decision maker can judge the worth of 
the objectives 7j(G). For an MOLQR formulation to be valid the designer 
must have enough understanding of each single objective in order to answer 
the trade-off question effectively. 

2.4 MOLQR formulations and examples 

The purpose of this section is to formulate some useful multiple objec- 
tive linear quadratic regulator (MOLQR) problems in which the designer 
can address the trade-off question. 


2.4.2 Integrated control system design 

Integrated control system design deals with the integrated design of 
subsystem controllers. The overall system is viewed as consisting of / sub- 
systems with the system state and control vectors z(f) and u(£) separated 
into subsystem state vectors x,(<) and control vectors u*(£), i = 1 
Each subsystem has a different objective function 

r oo 

MG) = / + uf (1)^(0] dt. ( 2 - 12 ) 

Jo 


2.4.1 The general linear quadratic regulator problem 
The integral quadratic objective function 

= r[x T (t)Qx(t) + u T (t)Rn(t)]dt (2-7) 

7o 

is a sum of terms 


± 7.; / M0-M0 ^ ij - 1, . 

Jo 

■ , n 

(2-8) 

±1 ’kt f Uk(t)m{t) dt fc t l= 1,. 

Jo 

. . , m 

(2-9) 


To solve this problem using the weighted sum method the single objective 
/ l 

J (G) = J2 u ’i J i(G), £> = 1, 0 < u/j < 1 (2-13) 

i=i »=i 

which is of the form (2-7), is minimized subject to any constraints on the 
control law. The weighting matrices Q and R take on the form 


Q = block diag{tiqQj, w 2 Q 2 , • • • , wiQt] 
R = block diag{uq R\ y w 2 R 2x ■ • 
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Fennel ami Black compute a single noninferior point for an integrated 
airframe and propulsion control problem using the weighted sum method 
[23). The weights were chosen to be the inverse of the optimums of the 
single objective problems, i.e., w, — 1/7* where 7* is the optimum for 
the single objective problcmi of minimizing 7,. The objective function wjls 
optimized for a full state feedback structure by solving a Riecati equa- 
tion. Suboptimal solutions of different structures, including a decentralized 
structure, were obtained by zeroing elements of the optimum full sLale feed- 
back gain matrix. Caglayan, llalyo, and Broussard reformulate Fennel and 
Black’s integrated airframe and propulsion control problem as a discrete 
time stochastic problem [ 10] . The weighted sum of objective functions is 
optimized subject to a decentralized control law using an iterative algo- 
rithm developed in [10]. However, neither of the above references ([10], 
[23]) investigate the trade-offs between airframe and engine performance, 
or the trade-offs between controllers of different, structures. This was done 
in [27] using the modified algorithm. The trade-offs are summarized in fig- 
ure T Significant cost, reduction was obtained using the modified algorithm 
compared to the suboptimal techniques in [23], table 1. 

2 1.3 Multiple regime conti itl system design 

in a mult iple regime control system design problem it is assumed that 
the plant may he any plant in a finite set of t time-invariant plants. The 
problem is to find the best control law that simultaneously optimizes the 
performance objectives of each plant in the set. Kadi plant has a different 
objective function 


■/,(<••')= / lzT(l)QiXi{l) + nl(t)Hinill)]dl (' 2 - 15 ) 

7o 

where Xi(t) and «,•(<) are the state and control vector for the i^ 1 plant. To 
solve this problem using the weighted sum method, vectors x(f), t /(<), and 

u(t) are defined which consist of the vectors Xj(f), i/,(t), and w»(t), of the / 
plants. The scalar objective function 7(G) then satisfies (2-13) and (2-7) 
with weighting matrices 


minimizing a model-following term subject to control energy con- 
straints on each of the helicopter’s two inputs. In the second example 
a control system for an aerodynamically unstable aircraft is designed by 
minimizing sensitivity subject to constraints on model-following errors and 
control energy. The trade-offs between model-following and sensitivity are 
examined 

2.4.5 Closed-loop trajectory sensitivity as au objective 
Consider systems of the form 


i(0 = A(p)x(t) + B(p)u(l) 

v(t) = < T(o 


(2-H>) 


where x(i), «(/.), and y{l) are defined earlier (1-1) and p 6 R r is a vector 
of parameters. To simplify notation only the single parameter case will be 
discussed below although the results easily extend to the multiple parameter 
case. ] h<* closed-loop trajectory sensitivity due to changes in a parameter 
p i may he reduced by adding the term 


dr. dr 

rr- QrTT-dt 
dp i <l}»\ 


( 2 - 2 ( 1 ) 


to the objective function (2-7) [20], [22]. The closed-loop trajectory sensi 
tivity function Aj (/.) = dx(t)/dpi is described by the state equations 


MO 

•Sj(0) = 0 


mp)Si(o + ™M x{l) + mp: 


dpi 


dpi 


At) + B(p) 


MV 

dpi 


(2-21) 


The term du(t)/8pi depends on the form of the control law. Fleming and 
Newrnann [20] augment the state vector x(t) with the sensitivity vector 
Sj(0- A full state feedback controller is then 

At) ~ -Frt*) “ *2Sj(0> (2-22) 


Q — lilock diag{uqQi , w 2 Qi, .... tviQi) 
R = block diagfu/j Ri , W 2 R 2 , ■ ■ - , • 


(2-16) 


in order to find the single best control law for all the plants in the set, the 
control laws of the plants must be constrained to be equal to each other by 
the constraint 


G = block diag{G ll G 2 ,...,G/} 
(j\ — G 2 — ' = G 1 . 


To implement this control law 5j (t) must be simulated. However when the 
control law (2-22) is substituted into (2-21) a term 0Si(t)/Dp\ appears. 
Fleming and Newrnann neglect this term. Hafez and Loparo [22] derive the 
necessary conditions for optimality for the control law 

At) = ~Fx(t). (2-23) 

Unfortunately the algorithm of section 1.3.1 cannot be directly applied to 
the necessary condition in [22] because the matrix Q is a function of the 
feedback gain matrix F . 


The plants in the set could represent different operating conditions or 
failure modes. The expected cost E{7} = could then be min- 

imized, where p\ is the probability of occurrence of plant i. In terms of 
decision theory, the problem is now one of decision making under uncer- 
tainty. The expected cost should be modified to include a weighting term 
w, for each term p,7,(G). The term Wi is determined based on the decision 
makers trade-offs. For example the decision maker may be unwilling to take 
the risk of incurring a large cost 7^(G) even if there is a small probability 
of occurrence of plant k. 

Fleming and Pashkevich use the goal attainment method to design a 
controller for a robot arm that operates in three different regimes [13]. The 
weighting coefficients Wi ~ 7/ , where 7,* is the single objective minimum, 
were chosen so that the same degree of under-attainment is achieved in each 
objective. In [27] a proportional plus integral controller is tuned for two 
possible plants, Gi(s) = (1/s + 1), G 2 (s) = 2/(4 s + 1). The objective is to 
minimize the model following error of the closed loop step input response 
of both plants, where the model is G m (s) = 36/(s 2 -f 7.2s 4- 36). A slowly 
decaying exponential is used to approximate a step input. The trade-off 
surface between model-following objectives 7!(G) and 7 2 (G) is shown in 
figure 5. The closed loop step responses for both plants for a weighting 
factor of .5 is shown in figure 6. Greater model-following for plant 2 may 
be obtained (at the expense of plant 1) by solving the problem for a weight 
of < .5. 

2.4.4 Model- following objectives 

If a model of the desired closed-loop plant is known, the model- 
following objective 


(Xrn(t) - x(0) T <2m(Zm(0 - x(t))dt, 


(2-18) 


may be used to minimize model-following error where x m (f) is the state 
vector of the model. The problem can be formulated by defining a state 
vector which includes both z(t) and x m (0- 


However, the problem can be put into a form in which the algorithm 
can be applied directly. This is done by augmenting the control vector ti(<) 
with du(t)/dpi and constraining the control law so that (2-21) holds. Given 
the control law (2-23) the formula for du(t)/dpi is 


du(£) 

dpi 




(2-24) 


Substituting (2-24) into (2-21) and evaluating at p = p° the state equations 
describing Si (7) become 


M0 = AS\(t) + Aix(t) - B\u{t) - Bu Sl { t) (2-25) 

where the input u 5l (t) = du(t)/dp x - -FS x {t) has been defined and the 
following notation has been used 


A = A(p°). 


At = 


9A(P) 

dpi 


lp° 


B — B(p°), 


, _ 9B(P) 

1 M P o 


(2-26) 


The problem is then formulated as an optimal output feedback problem 
with a constrained control law as follows: 


US) -u JKi&nn »)(.;%) 

*’ ='(*,< 0 ) 

( At) \-_(r oy 

V«M07 Vo FJKSiWj 


(2-27) 


Q- 




(2-28) 


The closed-loop trajectory sensitivity may be further reduced by in- 
cluding a dynamic compensator with state vector z(t) described by 


The ^-constraint method is used by Fleming [14] on two flight control 
system examples. A longitudinal controller for a helicopter is designed by 

ORKitN - 


(u(i)\_(E H\(x{t)\ 


OE poor yUAO, 


(2-29) 
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I li.-rrf* in- 


MV = jr Mt) f lf <Mt) 
<)}>\ #pi Oj) 1 


(2*ao) 


L lit* st ate vector must be further extended by the compensator stale trajec- 
tory sensitivity function S,(Q = () 2 (t)f<)p l , leading to the following stale 
equations 


•T (0 = 

rip I 


(2-ai) 


l lie slab? equation for JS' j ( / ) is found by defining the input 




y «(0 

<b\ 


ES x {l) + US tl {l). 


(2-:i2) 


The reformul.il ed problem for optimizing (2-7) plus (2-20) subject to I, lie 
control law (2-20) is: 
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Q = block fliag{Q. 0, Qi , 0} /? = block diag{ ft, ( /, <1, f /, } (2-34) 


It is important to note that in some cases greater insensitivity of the 
response due to a change in commanded variables might be obtained using 
upm-loop compensation as opposed to feedback control. The closed-loop 
trajectory sensitivity to additional parameters p, may be considered by 
extending the state vector with the sensitivity functions Si(t ). 

The closed-loop trajectory sensitivity to the relative position of the 
center of gravity of an unst able aircraft is investigated in references [21], 

[11], and [27], Results using the modified algorithm [27], for a full state 
feedback control law, are presented in tables 2 and 3, and figures 7 and 8. 
The objective function is J = (1 — w)J m + i uJ s + J u , where J m is the model 
following error, ./* is the sensitivity measure, and J u is the control energy. 
The main trade oir in the reduction of sensitivity is the increase in control 
energy. 

3.0 Summary 

In this paper the synthesis of a linear output feedback control law of a 
specified structure has been set in a multiple objective framework. The de- 
signer may trade-off any number of individual integral quadratic objectives 
such as state energy, control energy, model-following, trajectory sensitiv- 
ity, and subsystem objectives. Noninferior solutions may be generated by 
applying the modified descent Anderson-Moore algorithm of section 2.5 to 
a single objective problem that is the weighted sum of the objective func- 
lions. Dynamic compensators as well as static output feedback controllers 
can be designed using linear equality constraints on the parameters of the 
output, feedback matrix to impose the structural constraints. 
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TABLE 1 . - COMPARISON OF COSTS BETWEEN REFERENCE [23] AND 
CONSTRAINED OPTIMIZATION, FOR THE INTEGRATED AIRFRAME 
AND PROPULSION CONTROL EXAMPLE, a = 0.4848 


Control 1 er 
structure 

Reference [23] 

Constrained optimization 

■■ 

J a 

J e 

J* 

J a 

mm 

Centralized 

1367 

29.89 

2624 



— 

Decentral i zed 

3330 

29.40 

6422 


26.95 

2641 

st4( . 1 ) 

1759 

33.10 

3329 

■n 

42.37 

2763 


TABLE 2. - OPTIMAL COSTS FOR MODEL FOLLOWING VERSUS 
SENSITIVITY EXAMPLE OF SECTION 2.4.5 


w 

mm 



J 

0.0 

0.1 

5.2697 

.1277 

4.1116x1 0 — 1 ^ 
1 . 2278x 1 0“ 01 

7.767X10- 2 

2.086X10" 1 

4.1116x10-12 
1 .2303xl0- 01 


TABLE 3. - OPTIMAL GAINS FOR MODEL 
FOLLOWING VERSUS SENSITIVITY 
EXAMPLE OF SECTION 2.4.5 


w 

mam 


9*3 

B I 

O r- 

O O 

-1.4760 

-.0739 


-0.238 

-3.436 

-0.903 

-2.308 
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FIGURE 1. - OBJECTIVE FUNCTIONS J^g) AND J 2 <g) VERSUS 
SINGLE DECISION VARIABLE g. 



FIGURE 2. - NON INFER I OR SURFACE BETWEEN THE TWO OBJEC- 
TIVES J^g) AND J 2 (Q) OF FIGURE 1. 
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FIGURE 5. - NONINFERIOR SURFACE FOR NOMINAL VERSUS OFF 
NOMINAL PERFORMANCE EXAMPLE OF SECTION 2 . 1 . 3 . 
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(b) SENSITIVITY REDUCED DESIGN (w = 0.1). 

FIGURE 7. - TRANSIENT OF FLIGHT PATH ANGLE FOR DIFFERENT 
PARAMETER VALUES p. 
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(a) R ICC AT I DESIGN (w = 0). 
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(b) SENSITIVITY REDUCED DESIGN (w = 0.1). 


FIGURE 8. - TRANSIENT OF ELEVATOR DEFLECTION FOR DIFFER- 
ENT PARAMETER VALUES p. 
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